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ABSTRACT 

We study the — limit of the 0{N) Gross-Neveu model in 
the framework of the massless form-factor approach. This model 
is related to the continuum limit of the Ising model with random 
bonds via the replica method. We discuss how this method may be 
useful in calculating correlation functions of physical operators. The 
identification of non-perturbative fixed points of the 0{N) Gross- 
Neveu model is pursued by its mapping to a WZW model. 
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1 Introduction 



The aim of this paper is to discuss the critical regime of the two-dimensional 
random bond Ising model by applying non-perturbative methods to the mass- 
less phase of the 0{N) Gross- Neveu (GN) model in the limit iV — 0. It is well 
known (and briefly sketched below) that in the context of the replica method 
the analytic continuation ^ of the GN model describes the quenched 
averages of the Ising model in the presence of gaussian distributed random 
couplings . Perturbative calculations based on the GN model have 

then been extensively and successfully used for studying the behavior of cor- 
relation functions of the random model in the infrared regime where the GN 
model for < 2 is asymptotically free Being an asymptotically free 

infrared theory, the perturbative series is plagued though by the presence of 
Landau poles which do not permit to study the behavior of the correlators 
on their short distance scales and therefore to follow in general the change 
of the theory in passing from the short to the large distance scales. In this 
paper we show how in some cases we may get around this problem by using 
the non-perturbative massless Form- Factor (FF) approach proposed in [Q. For 
the two-point function of the energy operator of the random model we show 
for instance that the absence of Landau poles in the S'-matrix formulation al- 
lows us to have much more information on the correlator at intermediate (and 
small) scales than what can be obtained within a perturbative renormalization 
group approach. A more general theoretical problem consists in identifying the 
ultraviolet fixed point of the GN model in its entire massless range N < 2. 
At the end of this paper we present some considerations on a mapping of the 
GN model to an interacting Wess-Zumino-Witten (WZW) model which may 
be useful for further investigation of this problem. 



2 S'-matrix in replica space 



Let us begin our discussion with the continuum limit of the random-bond Ising 
model expressed in terms of a Majorana fermion (see for example for details) 
with the partition function given by 



Z[m{x)] = Dlijj] exp 



(fx ip — m{x)) ip 



(2.1) 



In the above formula, m{x) is a (position dependent) random mass term as- 
sociated to the lattice random bond interactions. We assume that m{x) has 
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a gaussian distributed probability P[m{x)] oc exp (— m^(x)/4(yf). Let us now 
make use of the standard replica method to compute the average of the free 
energy: 



Z 



N 



\nZ = lim . (2.2) 

N-^O N 

The quenched averages of correlation functions can then be obtained by adding 
a source term to the partition function and differentiating (|2.2|) with 

respect to it. The above procedure leads to an effective action described by the 
Gross-Neveu model with N ^ colors: 



S= (Fx 



N. (2.3) 



Irrespective of the number of colors N , the GN model presents at the classi- 
cal level an infinite number of local and non-local charges which are conserved 
just in virtue of the equations of motion 0,^. For N > 2, the integrability 
at the classical level is known to survive at the quantum level: the associated 
quantum field theory is asymptotically free in the ultraviolet region, massive 
[y] and presents in general quite a rich spectrum of bound states. The inte- 
grability of the quantum model implies the elasticity and factorization of its 
^-matrix amplitudes which can be then determined by employing the unitarity 
and crossing symmetry equations, together with the residue equations coming 
from the bootstrap principle ITUIP . 



For N < 2 the model is instead asymptotically free in the infrared region 
and presents at the quantum level a massless phase, as it can be argued (at least 
perturbatively) from the changing of the sign of the /3-function near g = ||12|| . 
Its spectrum consists in this case of only massless Majorana fermions which 
can be right- or left-movers with a dispersion relation parameterized in terms 
of a crossover scale M and the rapidity variable a.s E = p = (M/2) exp(6') anc 
E = —p = (M/2) exp(— ^) respectively [|1^]. Assuming that the integrability^ 
of the model also holds for A^ < 2, we propose the following exact ^-matrix 
involving its fermionic massless excitations (for the general formulation of mass- 
less scattering theories, see 

Sr,r = Sl^l = Sfi^L = '5' (2.4) 



^This statement can be actually proven along the arguments presented in if one is 
inclined to be rather cavalier about the meaning of the color indices in this range of N. 

^This essentially follows the original proposal made in |l3| for the massless S'-matrix of 
the GN model, with the difference that there the S'-matrix in the L — L and R — R channel 
was incorrectly assumed to be —1 on the basis of Feynman rules derived from the Lagrangian 
(|2.3|) which however involves for N < 2 irrelevant non-renormalizable interactions. 
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where 

S'^^i^e) = da^.d^^'a.^e) + 5:5ta,{e) + 5i5la,{e) (2.5) 

with 

m — 6 u 

rfi - ^)rfi + — irfi - ^ ^)r( — ^ + ^) 

^ r(^)r(i-^ + ^)r(i-^)r(i-^-^) ^ ' 

\2m' \2 2m ^ 2m' \2 2m' \ 2m 2ni ' 

and 

The above expressions formally coincide with those discussed in ||1T| for the 
massive phase of the model, the main distinctions being the different role played 
by the rapidity parameter as well as the different range of the parameter N in 
the two cases. The ansatz ( p.6|) is supported by the validity of the Yang-Baxter 
equations for all possible combinations of right- and left-moving particles as 
well as by the 0(iV)-invariance of the interaction. Notice that for N < 2 
the above amplitudes do not have poles in the physical sheet (a feature which 
is consistent with the massless phase of the model). For N ^ 2 they are 
continuous functions and at iV = 2 reduce to the S'-matrix amplitudes of the 
Sine-Gordon model at the Coleman transition point — > Stt 

(r2{e) = 0, 

rr (f)) - -2^^ r(i-^)r(i + ^) 

^^^^^ - ^ r(f)r(i-f) ' ^'-'^ 

^ 2m ' ^2 2m ' 



-2m r(i-/-)r(i + /-) 



\2mil V2 2mil 



3 Minimal Form-Factors 

Let us consider now the calculation of correlation functions by means of the 
spectral representation based on the FFs (see for instance [|6|, p!6| , p!7| for the 
relevant formulas and notation). For the two-particle FF the problem can be 
initially reduced to consider the three different channels whose S'-matrices are 
given by: 

CTiso = iVo"! + 0-2 + as , (T± = a2±0-3. (3.1) 
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For each of these channels, the minimal FF can be evaluated. As an example we 
construct the one for the isoscalar channel. Watson's equations for a massless 
flow are [HI (see also 



where a, = R,L and S stands for a scalar ^-matrix which in our case is 
(Jiso- In order to solve these equations it may be useful to employ an integral 
representation for aiso- 

aUe) ^ - exp I -2 / - smh ^- j | . (3.3) 

(Note that for = 1 this simply reduces to —1 - the S-matrix of the Ising 
model). Let us write Fa^ ^-Jy^^ as 

where gai,a2{^) satisfies 

gauaii^) = ~9a2,ai{-0) , gai,a2{0 + '^'^i) = 9a2,ai{-0) , (3.5) 

whereas fai,a2{(^) fulfills the equations 

. f 7dx e(^+i)--e- , fOx\] ^ , 
/„„„2W = exp I -2 y — -^^^ smh j | /„2,ai(-^) , 

/,„„2(^ + 27r«) = /«2,ai(-^^) • (3.6) 

The solution of (|3.5|) for the R — R and L — L channels is unique (up to a 
normalization constant) and is given by 

QrA^^) = 9LAd) = sinh (^^^ . (3.7) 

For the R — L channel any linear combination of e^ and e~^ is a solution. The 
exact computation of the two-particle FF for the energy operator for = 1 
(which corresponds to the pure Ising model) and its expected infrared behavior 
at A^ = fixes the solution in a unique way 

9rA0) = • (3-8) 
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As for ( |3.6| ) we take the same solution in the R — R, L — L and R — L channels 
(this assumes that fR^hiO) = fL,R{0)): 

fR,R = fL,L = fR,L = f , (3.9) 
with the function / given by 



f{e) = exp -2 / — 

J X e^ + 1 smh x 

K 



(3.10) 



This can be also expressed as 



oo 



^^^^-n(r(2 + .)r(i-A + ,)j 

r(|-^ + fc)r(f + ^ + fc)r(i-A__^ + fc)r(-A + -^ + fc) 
r(2-^ + fc)r(i + ^ + fc)r(|-A__i_ + A;)r(i-A + _^ + fc)- 

Using its infinite-product representation one can show that f{6) satisfies the 
following functional relation: 

which may be useful in the computation of higher-particle FFs. The constant 
C\ is given by 

^ r(|) - / r(l- A + A:)r(| + fc) y (l + A:)(| + fc) 

' r(-A)r(i - A) n l^r(2 + A;)r(i -i^ + k)) (-a + fc)(i - a + fc) ' 

(3.13) 

its numerical value at = is = 0.5854 . . .. 



4 Two— point function of the energy operator 

Let us specialise our analysis to evaluate the (average) correlation function of 
the energy operator (e(x)e(O)). In terms of the replica, this is equivalent to 
evaluate the correlation function 

1 ^ 

T^E (efeWe;^(0))U=o (4-1) 

k=l 
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in the GN model. At the infrared fixed point, the energy operator is given in 
terms of the fermions by ek{x) ~ —iiplr\x)iljl^\x) where iIjI^\x) and iIjI^\x) 
are the chiral components of the original Majorana fermions. Duality and spin 
reversal symmetry translate into invariance of the Lagrangian (|2.3|) under the 



two transformations tp^^^ t-^ "^i^"* ^ind tp^ ^ t-^ ipl or tp^^^ ^ "^k^^ 

^ ^ —ipl. \ Under these two transformations, the energy operator changes 
its sign. Since these discrete symmetries are preserved by the perturbation, one 
concludes that the only non- vanishing FFs of the energy-operator are those with 
an odd number of both left- and right-moving particles. In particular, the first 
non-trivial FF is 

f5\9„92) = (0| (E e^(0)) m9,)L,{92)) . (4.2) 

Here \Ri{9i)Lj{92)) corresponds to an asymptotic state of one right-moving 
and one left-moving fermionic particle. 0{N) and Lorentz invariance fixes this 
FF (up to an overall normalization constant) to be given by F^J^{9i,92) = 
Si,jFii^L{9i - 6*2) with Fji^L as in eqs. (Q, (p?8|), ( p.lO| ). In the previous ex- 
pression there could be in principle a supplementary factor (a 2m periodic 
symmetric function without poles, i.e. a polynomial in e^'). This factor is ac- 
tually absent in the pure Ising model (A^ = 1) and its absence will be further 
justified later on by the correct asymptotic behavior of the correlation function 



The contribution of ([4.2|) to the two-point correlation function of the 
energy operator is given by 



(e(x)e(0))(2) ~ J ^ J ^ \Fn,L{9i - 92)\' exp ' 



= Cf\Mr) (4.3) 

where the superscript indicates that this is the two-particle contribution and 
X = [ir, 0) i. With the change of variables 'j = 9i + 92, 9 = 9i — 92 we obtain 



00 , „ 00 



Cf)(Mr) = 1 I ^ I g \FrA9)\' exp (-Mr ei cosh £ 

--00 —CXD 

= [ KoiMr el) \Fn,Li9)\' ■ (4.4) 



''Actually, there is another two-particle contribution to this correlation function, namely 
the one coming from Fl^r. However, it gives the same contribution as Fh^l. 
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Using ( |3.4| ), ( |3.8| )-( ]3rT0| ) we can express |-Fr,l(6')P as (from now on we will 
specialize our formulas to the limit N ^ 0) 



R,L 



1 + 4 



<l>(^) (4.5) 



with 



<l>(^^) = exp |2 7^e-^ tanh^ cos^ | . (4.6) 

It is straightforward to evaluate this last integral numerically, it is bounded 
from below by and from above by 



oo 



2 J ^e-" tanh^ (^0 = 0.315384 . . . , (4.7) 



i.e. 

1 < $(0) < 1.370786 (4.8) 

for any value of 6. In particular, the factor $(^) will not influence the 6 — >■ 
±oo asymptotics of the FF. Thus, we can replace ^{0) by a constant in order 

to discuss the infrared behavior of the contribution ( [4.4|) to the correlation 
function. Putting 

$(e) = l, (4.9) 

we find 



Cp)(Mr) =7^ , ^ , Ko{Mr p) , Mr > 1 (4.10) 



dp 

where we have made the change of variables p = . Upon expansion of the 
integrand, the leading infrared behavior is found to be given by 

^"''''"'•>- 2.(Afr?^lnAfr (^ + "(n7I?F)) ' '"^) 

This fits nicely with renormalization group results, as the following discussion 
shows. At the one- loop order, we have (see e.g. eq. (4.10) of ref. 0) 

{e{x)em = r-' + 0{g,) (4.12) 

1 + 87r5'o In f 
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0.01 



Cp\Mr) 



0.0001 
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Figure 1: Two-particle approximation to the correlation function of the energy 
operator (full line), the one-loop perturbative result (long dashes) and the 
two- loop result (short dashes). 



and including the two-loop contribution one arrives at 

4^0 



{e{x)e{0)) = r" 



(l + 87r^o lnr)(l-27r^o 
4vr5fo 



X 



(4.13) 



1 + 



1 + Stcqq In r 



[In (1 + Svr^olnr) - 47r^olnr] + 0(^^ 



where Qq is a bare coupling constant. Here we have chosen the normalization 
such that both (|0^ ) and (|]T3|) agree with { ^AA\ ) at M = 1 for r oo. Setting 
M = 1 we obtain a good agreement with ( [4.4|) in the region Mr > 10 for ( [4.12 ) 
with qq = 0.23 and for ( [4. 13] ) with qq = 0.038. One can see that the Landau 
pole in the two-loop formula is shifted towards smaller scales and the behavior 
of (|4.13| ) is more regular around Mr = 1 in Fig. 1. Higher loop corrections 
in the pertubative framework are expected to move this pole to smaller and 
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smaller scales. The main point is the absence of this kind of singularity in our 
result based on the FF approach. 

5 Comments on higher-particle FFs 

Apart from logarithmic corrections, the energy operator has canonical scaling 
dimensions in the infrared and, as explicitly checked for similar operators in 
other massless scattering theories, a rapidly convergent spectral series for its 
correlation functions is expected in this case P,p!9| . This expectation seems in- 
deed confirmed by the above calculation. Under these favorable circumstances, 
one is essentially released from computing the more complicated and cumber- 
some expressions of the higher-particle FFs for most reasonable purposes. It 
should be pointed out though, that these favorable cases do not always occur. 
Consider, for instance the (averaged) two-point correlation function of the spin 
operator. In this case two kinds of technical problems arise: the first difficulty 
is that the integrals entering the spectral series ought to be regularised in or- 
der to cure the infrared divergencies. This is however a well-known problem 
of the spectral representation based on massless FFs , a problem that now 
can be easily handled following the suggestion of ref. |2^. The second diffi- 
culty is the most serious one, namely the apparent necessity of employing in 
this case all higher-particle FFs of the spin operator in order to recover both 
its infrared and ultraviolet behavior. The problem is particularly difficult here 
since there is a huge difference in computing higher-particle FFs in diagonal 
rather than in non-diagonal scattering theories (in the diagonal case one is of- 
ten able to obtain a closed expression for all the higher-particle FFs, see for 
instance ^^). Although progress has been recently achieved on the FF problem 
in some non-diagonal scattering theories |T8| , p2| , the general determination of 
the higher-particle FFs for non-diagonal S'-matrix models still remains an open 
problem of a (quite) mathematical difficulty. It would be clearly interesting to 
try to develop the massless FF approach further in order to deal successfully 
with this class of operators. 

6 Non-trivial fixed points 

In the last part of this paper, we present some considerations which may be 
useful to better understand the structure of the fixed points in the 0{N) GN 
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model and in particular to determine its ultraviolet fixed point in the range 
N < 2. The theory behaves once again differently for N > 2 and for N < 2. 
More is known about the structure of the fixed points of the former case rather 
than of the latter. Let us in fact recall that for N > 2 the model is massive but 
ultraviolet asymptotically free. This means that at short distances it reduces 
to free Majorana fields (with central charge Cuv = N/2) whereas at large 
distance scales no massless degrees of freedom are left and its central charge 
is therefore Qr = 0. Conversely, for N < 2 the model is asymptotically free 
in its infrared scales and there we have correspondingly Qr = N/2. However, 
since in this case the model is massless along all its flow from large to short 
distance scales, the central charge in the ultraviolet limit is one of the dynamical 
data which remain to be determined. Some insight can be gained by rewriting 
the partition function of the model defined by the action (|2.3|) in a way more 
appropriate for studying the existence of some other fixed point in addition to 
the (trivial) one at g = {a similar procedure has been used in to study 
non-trivial fixed points in the chiral SU{N) GN model). To this end, we first 
make use of the identity 

(V'a^a)' = -2(V'a7M^>b)', (6-1) 

where are the 0{N) generators normalized as trT^T^ = 6^^ (this identity 
is a direct consequence of the fact that we are working with Majorana fermions). 
Then the quartic term in the interaction can be traded for a quadratic term 
via the introduction of an auxiliary field through the identity 



exp \^2g J (fx (jJal^Jab'^}^ 

I DA^ exp (^^ / dMA^r + 1 1 d'x{^P^,TU)A^^ , (6.2) 
and hence the partition function can be written as 
Z = J Dij^DiJaDA^exp - J fx Ua<^6ab + i 4aMh + ^K^i 



. (6.3) 



The key observation now is that the fields in eq. ( |6.3|) can be decoupled through 
the transformations 

A = -dhh^, A = -dgg^, ^ = {l (6.4) 
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where h and g are 0{N) matrix- valued fields, and ^ stands for transpose. 
Taking into account the Jacobians of the above transformations [2^] (for more 



details compare also Section 9 of [^), we can rewrite the whole partition 
function product of decoupled sectors 

Z = ZfiZgi^Zi^t , (6.5) 

where Zg is the partition function for the free Majorana fermions, the 
ghost partition function arising in the computation of the Jacobians and 

Z,nt = J DhDgexp{{l + 2C,) {T[h] + r[/]) 

- (^i^|^ + ^)/rf^xtr(9/./.W)}, (6.6) 

where a is a parameter that keeps track of regularization ambiguities, C^j is 
the dual Coxeter number of 0{N) and r[u] is the WZW action |2^. Hence, 



modulo decoupled conformally invariant sectors, we have an effective theory of 
interacting WZW fields. 

One comment is in order regarding the regularization ambiguities arising in 
the evaluation of the Jacobians. The quadratic term in in ( p.3|) breaks both 
gauge invariance and local chiral invariance, hence there is no a priori reason to 
choose a regularization preserving any of these two symmetry transformations. 
We will fix the parameter a later on. Using the Polyakov-Wiegmann identity 



[ p^ it is straightforward to see the existence of a fixed point given by the value 
of the coupling constant gi = 7r/(4(l + 2Ct,)(l — a)). The effective partition 
function at this value may be written as 

Zi„t|,, = / D~gexp ((1 + 2C,)r[g]) , (6.7) 

where the identification g = g^h has been made and the integral over h has 
been factored out. Taking into account the free fermions and the ghosts, the 
partition function ( |6.5[ ) corresponds to a conformal field theory whose Virasoro 
central charge is given by 



ci = — - N{N + 



{1 + 2C,)N{N -1) 



N / iV- 1 , 



2(1 + C'„) 

There may be however another fixed point, given by the value g2 = — 7r/(4(l + 
2C„)a), where we have 

2i„tU = J DhDgexp{{l + 2C,)(T[g^] + T[h])}, (6.9) 
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and the corresponding Virasoro central charge is given by 



C2 = — -iV(iV-l) + 2 



[1 + 2C,)N{N-1] 
2(1 + CJ 



We have now to choose the parameter a and secondly to identify the fixed 
points in the two ranges N > 2 and N < 2. 

Let us first see how the above points can be settled in the well-understood 
case N > 2. Since in this region the model is asymptotically free in the ultra- 
violet and massive otherwise, no non-trivial fixed point is expected. Therefore 
one should choose = cx) in ( |6.3| ) since then the integration over the gauge 
fields leads to constraints that eliminates all the degrees of freedom yielding a 
c = theory. This suggests the choice a = 1 {i.e. the gauge invariant regular- 
ization) for which we have an infrared fixed point in the strong coupling limit 
(?i — > oo. Plugging Cy = N — 2 in the corresponding expression ( |6.8D of the 
central charge, we find indeed the correct value ci = 0. Notice that with this 
choice of a, the value of g2 is negative and does not correspond to a physical 
fixed point, because this is incompatible with the unitarity of the GN model 
for > 2. 

The GN model at = 2 is somewhat special since it is equivalent to the 
massless conformally invariant Thirring model with c = 1 for —^<g<C!C. At 
N = 2 one has = and retaining a = 1 as before, one finds that (71 = 00 
still corresponds to a theory with ci = which is consistent with approaching 
N = 2 from above out of a massive regime. The other fixed point is now located 
at the end of the line where c = 1, i.e. g2 = —- and C2 = 1 which seems to be 
more appropriate for the approach to A^ = 2 from below out of the massless 
regime. 

Let us now discuss the fixed points in the regime N < 2 assuming that both 
formulas a = 1 and = N—2 also apply here. The GN model is generally non- 
unitary in this regime, a fact which opens the possibility to consider also a fixed 
point with a negative value of the coupling constant. With the above choice of 
a and Cy, the central charge for the fixed point gi is identically zero also for 
N < 2. Let us consider the second critical value of the coupling constant, i.e. 
g2. Note that g2 is negative for A^ > | but positive otherwise and therefore also 
compatible with the statistical interpretation as the random-bond Ising model 
for A^ — > 0. The corresponding value of the ultraviolet central charge is then 
given by ( |6.10D , i.e. upon inserting Cy = N — 2 

N(2N - 3) , , 

Cn.= ^^ -. (6.11) 
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The above formula does not seem to apply to the case = 1: in fact it predicts 
Cuv = — 1/2, but instead we expect to find in this 1/2 for the 

simple reason that it is impossible to construct a quartic fermionic interaction 
with only one Majorana fermion. The reason of this mismatch seems somehow 
interesting: first notice that using (|6.1CI|) the variation of the central charge 



from the short to the large distances reads 

Ac = iV(iV-l)^^ . (6.12) 

It is a general result of two-dimensional quantum field theories that such vari- 
ations of central charges satisfy the sum rule |2^ 



Ac = ^ J d^x\x\^ (e(a;)e(0)) , (6.13) 

where G(x) is the trace of the stress-energy tensor. This operator is propor- 
tional to the quartic interaction term in the Lagrangian ( p.3|) (the propor- 
tionality constant being the j3{g) function of the model). Therefore, the term 



A^(A^ — 1) in (|6.12|) has a pure combinatorial origin. This factor itself of course 
vanishes for = 1. However, after pulling out this combinatorial term, what 
is left in (|6.12|) or in (|6.13|) may be interpreted as "Ac per unit of replica", i.e. 



a quantity which still depends on but which has lost any reference to the 
color indices of the theory. Let us denote it by D{N). For = 1 this quantity 
can be easily computed by means of the Green function of the free fermion 
resulting in a (logarithmically) divergent expression 

D{1) J d'^x \x\^ . (6.14) 

For A^ 7^ 1 we expect that the divergence is cured bv the presence of interactions 
so that we expect D{N) to be a finite quantitya. The above considerations 
suggest then for A^ = 1 there may be a sort of anomaly which conspires to 
give a non-zero value for Ac. In this case the correct procedure might be to set 
A^ = 1 in ( |6.10|) before inserting the expression of Cy, i.e. to use a regularization 
expression for D{N = 1). 



^At iV = 2 it is the vanishing of the f3 function which is responsable for the vanishing of 
Di2). 
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7 Conclusions 



In summary, in this paper we have apphed massless FFs to the critical regime of 
the random-bond Ising modeL Compared with the perturbative approach, the 
energy correlation function obtained with this method is well behaved in the 
whole range of scales. The technical difficulties encountered in the computation 
of higher-particle FFs can possibly be circumvented by using methods similar 
to the ones used in |18,2^, where FFs for some non-diagonal S'-matrices have 
been computed. One motivation of our approach is to open the way for a more 
deep understanding of the ultraviolet behavior of the massless flow of the Ising 
model from its pure to its disordered regime and, more generally, for all the 
0{N) GN models with N < 2 which appear in the replica approach of the 
original random problem. To this end we have also pointed out a possible way 
of reaching the ultraviolet fixed point by a mapping of the GN action to a 
WZW model, relating the strong and weak coupling regimes. This approach 
predicts the presence of a non-trivial fixed point with central charge given for 
7^ 1 by ( |6.11|) . The lacking of a sound mathematical definition for the 
relevant quantities of the group 0{N) for N < 2 makes it highly interesting to 
have independent information on this issue. In this respect, the most natural 
approaches are those based on the Thermodynamical Bethe Ansatz or the 
afore- mentioned c-theorem sum rule p8[. The application of the two approaches 
seems however presently obstructed by technical difficulties related both to the 
non-diagonal nature of the S'-matrix and to the subtle problem of how to take 
the analytic continuation of the mathematical expressions for continuous values 
of N in the range N < 2 (this is particularly severe for a Thermodynamical 
Bethe Ansatz approach). The solution of these problems together with the 
massless FF approach proposed here may give interesting non-perturbative 
information in the field of disordered systems. 
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